A Density Matrix Renormalization Group study 
of Excitons in Dendrimers 
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We introduce the density matrix renormalization group (DMRG) method as an efficient com- 
putational tool for one-exciton approximations with off-diagonal disorder. This method allows 
us to reduce the computational effort by targetting only a few low-lying eigenstates at each 
statistical samplings, in contrast to the exact diagonalization methods that compute the whole 
spectrum. As an application of the method, we study excitons in two families of branched 
molecules called dendrimers using a recently introduced simple model. We compute the ab- 
sortion peaks for these dendrimers varying their generation number g and number of wedges 
w. 

PACS numbers: 75.10.-b, 05.50.+q, 71.35. Cc 



I. INTRODUCTION 

Dendrimers are large, highly branched polymers |lj] 
(see Fig.Q). Some individual dendrimer molecules ex- 
ceed ten nanometers in diameter. These hyperbranched 
molecules are composed of a central core to which repet- 
itive dendritic branches are attached (2j. Dendrimers 
have several distinctive architectural components. Some 
of them are: i) an initial core also called focal point; ii) 
interior layers called generations, composed of repeating 
units radially attached one after another. This number 
determines the branch length; hi) the degree of connec- 
tivity of each single molecule in a given site of the macro- 
molecule, i.e., the number of nearest-neighbours for a 
given site in the associated lattice model. This can vary 
from site to site; iv) the number of initial strands at- 
tached to the focal point can be made variable and leads 
to various separated branches that are called wedges or 
lobes. In other words, this is the core multiplicity. Thus, 
we can partially characterize the lattice model of a den- 
drimer using the notation 



B(w,c,g;N) 



(1) 



to denote a dendrimer with w wedges, connectivity c, g 
generations and total number of sites N . The number of 
sites can be obtained from the previous data (w, g, c). In 
the following we shall use this nomenclature and a planar 
representation of dendrimers (see Figs.|[ |||). Particular 
examples of dendrimer lattices are a binary tree, a Bethe 
lattice or a Cayley tree. 

The field of highly branched macromolecules can also 
be addressed from the viewpoint of polymer chemistry. 
We may classify polymers in two big groups: linear and 
branched polymers. Standard polymers elongate in a lin- 
ear fashion §,§,§], while the new dendritic polymers 



are branched M, Q, ||. Dendritic polymers differ from 
linear polymers in that the latter consists of long chains 
of molecules, like coils, crisscrossing each other, while 
dendrimers have numerous chain-ends that can be func- 
tionalized. Because of this, dendritic molecules can be 
constructed with discrete domains having different prop- 
erties. With a quite descriptive and friendly terminology, 
standard polymers are called "spaguetti-like" while den- 
drimers are called "meat-balls" . 

The hystorical origins of dendrimers are in the quest 
for large, substrate-selective ligands, in which several 
research groups became interested yielding the synthe- 
sis of "tentacle" || and "octopus" jl0|, |ll| molecular 
compounds, where long branches radiate from a cen- 
tral hub or a macrocycle. The first who reported a 
dendrimer synthesis was Fritz Vogtle in 1978 JlH , in- 
troducing the concept of cascade reaction and cascade 
molecules, later known as dendrimers. Then, research 
groups led by Tomalia |l3| and Denkewalter [jl4| de- 
vised routes whereby stepwise polymerisations could be 




FIG. 1: A pictorical representation of a dendrimer with 
parameters w = 3, c = 3, g = 7. 
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FIG. 2: Schematical representation of the wedges in a 
dendrimcr with a core with 6 wedges. The wedges are 
ascmblcd together as in the convergent method of sinthe- 
sis [[15). 

achieved, producing highly branched polymers with ex- 
tremely low polydispersities. 

Dendrimers are constructed from branching units and 
a core. Their synthesis begins by attaching branched 
molecules to an initial core structure. Iterative poly- 
merization of more branched molecules ultimately leads 
to a globular structure unable to accommodate further 
branching because of steric hindrance. Using multistep 
repetitive syntheses, chemists can construct the den- 
drimer layer by layer (generation by generation). This 
method is known as the "divergent method." It was de- 
veloped in 1978 by Fritz Vogtle. It is a difficult and 
tedious method for their synthesis, but it does allow 
exquisite tailoring of architectural features. Later, other 
alternative methods were proposed such as the "conver- 
gent method" , were dendrimers are built from the sur- 
face to the core, again using a reiterative sequence of 
reactions. This method was first introduced in Jl5| . 

Dendrimers may have a huge variety of constituent sin- 
gle molecular groups. Typically they contain carbon- 
based organic molecules such as polyamidoamines, 
aminoacids, DNA, and sugars, but also they are being 
constructed from organosilicons and organic/inorganic 
hybrids. 

There are many examples of dendrimers whose consti- 
tutions have been designed with a purpose in mind. By 
an accurate choice of the building blocks to be used and 
specific functionalization in the periphery (external sur- 
face) , many of their properties can be controlled such as 
the molecular weight, size, shape, chirality, density, vis- 
cosity, polarity, solubility, flexibility and surface chem- 
istry of the resulting macromolecules. This leads to new 
materials with a great potential for future applications. 

Interest in dendrimers has mushroomed after produc- 
ing the wide range of structures that are known today, 



and a great deal of potential applications for dendrimers 
have been proposed ||. Among the many current 

and future applications we may cite drug delivery sys- 
tems (in chemotherapy for instance), gene therapy, ma- 
terials science, magnetic resonance imaging (MRI) con- 
trast agents, molecular recognition, antiviral treatment, 
energy converters, etc. In material science, dendrimers 
may themselves act as building blocks in the fashioning 
of larger supramolecular structures. 

Dendrimers include a great volume in relation to their 
molecular weight as a consequence of their dendritic 
structure. They may quite well be compared to bushes 
and trees in nature (see Fig.|l]). Cavities inside them can 
host guest species such as ions or another molecules. The 
potential of dendrimers in the fields of host-guest chem- 
istry and nanotechnology relies on the subtle engineering 
of their architectures to pre-defined designs and taking 
advantage of these cavities inside them to encapsulate 
the guests. 

Of all the many applications of dendrimers, we shall 
focus on their optical properties that confer them the 
potential applications as light-harvesting molecules, i.e., 
molecular antennas. There have been recent experiments 
addressing these issues by studying electron transfer in 
dendrimers of several types and their optical properties 
for absortion and emission of light Jl6| . 




FIG. 3: An example of extended dendrimer with gener- 
ation number g = 4 and wedges w = 3. At the ramifi- 
cation sites the coordination is c = 3 like in the compact 
dendrimers, while the new feature is the presence of in- 
termediate sites at the legs with c = 2. 

The functioning of dendrimers as single molecule pho- 
tonic antenna systems relies on the large number of chain- 
ends that can be functionalized to absorbing photons at 
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the periphery and a subsequent efficient transfer of the 
absorbed energy to the center where a fluorescent trap 
or a laser dye can be placed at the dendritic core 

Following these experimental studies, we shall focus for 
the moment in two types of dcndrimcr familcs: compact 
dendrimers and extended dendrimers. They are shown 
in Fig.|^. These are phenylacetylene dendrimers where at 
each node or site of the dendrimer lattice there are ben- 
zene molecules. Depending on the position of these ben- 
zenes, the branching of the site can be either two (para 
position) or three (ortho position). Compact dendrimers 
consists of three Cayley trees joined together at the focal 
point. This implies that the number of wedges is w = 3. 
All their sites exhibit ortho positions leading to a connec- 
tivity c = 3. The number of generations is variable. In 
p6| dendrimers up to generation 5 were studied. Thus, 
a compact dendrimer is a D(3,3,5;iV(#))-dendrimer (@) 
with 

g 

N(g) = 1 + 3^2"^ = 3 x 2 9 - 2. (2) 

71=1 

This is an exponential growth. The first elements of this 
series are D4, D10, D22, D46, D94, etc. The experimen- 
tal studies of pl| show that the electronic excitations of 
these fractal structures are of a localized nature. 

Extended dendrimers have also a threefold symme- 
try around the focal point, but they have generation- 
dependent segments lengths. Namely, their sites can 
be of variable connectivity, cither c = 2 (para position) 
or c = 3 (ortho position) depending on the generation. 
Thus, an extended dendrimer is a D(3, [2, 3], g; N(g))- 
dendrimer (til) with 



Compact Dendrimers < — ► Gapped Spectrum 
Extended Dendrimers < — > Gapless Spectrum 

A complete ah initio computation of the electronic 
properties in these two families of dendrimers would re- 
quire to include many different effects and a truly many- 
body model for these supermolecules. However, recently 
a very simple phenomenological model of one-exciton 
processes has been proposed by Harigaya ^||, |2i[| 
which produces good fits for the experimental absortion 
peaks in extended dendrimers. 

In this paper our aim is to propose the density matrix 
renormalization group (DMRG) method as a good com- 
putational tool for dealing with one-exciton process in the 
presence of disorder. This method allows us to perform 
an exhaustive study of excitons in the Harigaya model to 
be explained below, for both families of dendrimers. 

Our studies are theoretical, we consider dendrimers as 
a useful model for investigating the dependence of phys- 
ical properties on molecular size and topology. This pa- 
per is organized as follows: in Sect. II we introduce a 
simple Frenkel Hamiltonian describing exciton processes 
and concentrate in one-exciton approximations using a 
recently introduced model Hamiltonian; in Sect. Ill wc 
propose the DMRG method as an appropriate computa- 
tion method for excitons with disorder in extended and 
compact dendrimers, and we present our numerical re- 
sults in a variety of situations with different numbers of 
generations g and wedges w. Sect. IV is devoted to con- 
clusions and future prospects. 

II. EXCITONS IN DENDRIMERS 



N(g) = l + 3^(.9-n)2 n - 1 +3 2^ 1 

n=l 

= 9 x 2 9 - 1 - 3.g-2, 



(3) 



where the factor in (|3|) l(g, n) = g — n+ (n = 1, . . . , g) is 
the length of each segment in the dendrimer. This is also 
an exponential growth. The first elements of this series 
are D4, D10, D25, D58, D127, etc. 

Despite this small difference in the structure of ex- 
tended dendrimers, it has important physical conse- 
quences which have been reported experimentally [ff6| : 
electronic excitations become delocalized with increasing 
size (generation) of the supermolecule. This experimen- 
tal evidence comes from the study of dendrimer absortion 
spectra, which clearly shows that the first absortion peak 
for the compact family remains constant at a given wave 
length, while in the extended family this wave length dc- 
cresase with increasing size. In other words, in extended 
dendrimers the gap for excitations from the ground state 
to the first excited states tends to close while it remains 
constant in compact dendrimers. Thus, there is experi- 
mental evidence for the following scenario: 



A consequence of the aforementioned experimental 
studies [[l6) is that the family of extended dendrimers 
can serve as artificial light-harvesting antennas. This 
has been demonstrated experimentally jL8j. The light 
is harvested over a wide area on the surface of the den- 
drimer and funncled by dipolc-dipole interactions to the 
single active site at the focal point where energy conver- 
sion takes place. Based on the optical spectra of the first 
absortion peak |16 , Kopelman et al. conjetured that 
the electronic excitations in extended dendrimers are lo- 
calized on the linear segments of the branched molecules. 
Theoretical studies have confirmed this fact showing that 
the relative motion of photogencrated electron-hole pairs 
(excitons) is confined to the various segments and energy- 
transfer may then be described by the Frenkel exciton 
model [Hi). 

In the Frenkel exciton model, at each site of a lattice 
there is an active center modeled with a two-level system. 
At half-filling, there is an electron per site on average 
and in the absence of interactions among centers, the 
ground state |GS) is formed by putting each electron at 
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FIG. 4: Simplified treatment of the phcnil rings as single 
sites. 

the corresponding ground states |vac) of the two-level 
systems, that is, 

|GS) = | vac) 1 1 vac) 2 . . . | vac) at (5) 

A local excitation is formed by promoting one of the 
electrons at a given site to the first excited level |exc). 
This is a localized electron-hole pair that models an ex- 
citon. An example of this is the state 

| Ex) = |i) = |vac)i . . . |exc) x . . . |vac)jv- (6) 

Let us introduce fcrmionic creation and annihilation 
operators c\ a , Ci^ a for these localized centers, where i de- 
notes the lattice site and a is a label that stands for the 
energy levels: a = for |vac) and a — 1 for |exc). Then, 
the localized energy levels can be represented as 

| vac), = ct |0) 
|exc)i = c{ JO) 

where |0) is the vacuum of the fcrmionic operators. The 
excitonic operator that creates one exciton from the lo- 
calized ground state |vac)j is denoted by b\ and acts as 
follows 



where E{ are chemical potentials that represent the en- 
ergy gaps at each site of the lattice A and Jij is a hop- 
ping integral for nearest-neighbour centers that in- 
troduce interactions between the site centers [gl]. The 
origin of these couplings is the residual electric dipole in- 
teractions between the molecules modeled by the active 
centers. It is assumed that this interactions are short- 
ranged. Their effect is to allow excitation transitions 
between nearest-neighbours sites. 

Harigaya has recently proposed to study one-exciton 
processes in the family of extended dendrimers introduc- 
ing off-diagonal disorder |^2|, |2^, [24|] . This simple model 
produces a good accurate phenomenological fit to the 
experimental data of the first absortion peak reported 
in Jlijj . Harigaya considers a first simplification of the 
molecular estructure in extended dendrimers in order to 
make the analysis more tractable. This simplification 
amounts to replace the complicated phcnil rings by sim- 
ple activation centers as depicted in Fig.|[ A second 
aproximation considered in |^2|, |2^, ^4) is to restrict the 
study to the one-exciton aproximation within the Frenkel 
Hamiltonian (|^) for excitons. We shall refer to this sim- 
ple model based on these two aproximations as the Hari- 
gaya model. A further assumption in this model is the 
presence of off-diagonal disorder. 

To derive the one-exciton approximation of the Frenkel 
Hamiltonian (^|) we project this Hamiltonian onto the 
Hilbert space of one-excited states (g) by considering a 
general superposition of localized excited states as follows 

\^) = J2 c -\ x ) ( 10 ) 

where C x are unknown amplitudes of the one-exciton 
states. Next, we set up the Schrodinger equation in this 
sector 

# l^cx) = £cx|*cx) (11) 

The solution of this equation is equivalent to solv- 
ing the Schrodinger equation for an associated one-body 
Hamiltonian Hi cx that takes the following form 



|exc)i = &]|vac)i, with b\ = ct^o- (8) 

The excitonic operator b\ creates a localized electron-hole 
pair at the site i and has bosonic character. 

When interactions among centers are brought to the 
system, the interesting question is to see how one these 
excitons evolve under quantum fluctuations. A simple 
model for these interactions is presented in the Frenkel 
Hamiltonian, namely, 

H = Y^Ei b\b t + £ J lt j (b\b 3 + b)h) (9) 



ffiex = J2 E *M\ + + h - c -) (12) 

where \i) are one-particle states representing an exciton 
state at site i. That is, an onsite exciton state is assigned 
to each phenyl ring (see Fig|| and Fig.|J). 

In extended dendrimers, interaction strengths between 
neighouring dipole moments may vary among position of 
dipole pairs. In the Harigaya model it is assumed that 
they are randomly distributed and the following Gaussian 
distribution function is chosen for the coupling strenghts 

P(J,,) = ^_ e -^) 2 (13) 
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where J is the standard deviation of the interaction and 
the mean value of the interaction is taken as zero. 
Further, the excitation energies Ei are also assumed to 
be uniformly distributed 

Ei = E, Vi (14) 

Thus, the two adjustable parameters of the model are 
E and J. The former specifies the central energy position 
of excitons in the optical spectra. According to |^2|, |2^, 
pif , we shall take the values of these parameters as E = 
37200cm -1 and J — 3552cm -1 , which were found to give 
the best fit to the experiments |l6| ]. 

The diagonalization of ( |l2| ) under these assumptions 
gives the energies of one exction states measured from 
the ground state. The energy position of the optical ab- 
sortion edge is always given by the lowest eigenvalue be- 
cause the state with the lowest energy is always allowed 
for dipole transition from the ground state ^3], p4[ . 
Thus, the values of the absordion edge E a b (energy of 
the lowest optical excitation) for the family of extended 
dedrimers is computed as 



E, 



ab 



E-(e )J 



(15) 



where (eo) is the average of the ground state energies 
eo obtained by diagonalization and sampling of the one- 
exciton Hamiltonian in eq. (|l2]). 

To solve this model of stochastic fluctuations in the 
Hamiltonian matrix elements that simulates the random 
distribution of couplings, in |25 
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an exact diag- 
onalization procedure was employed for a sampling of 
10000 realizations. This was carried out up to the fifth 
generation of extended dendrimers which has D127 sites. 

As in the computation of the first peak of the absortion 
spectra only the computation of the ground state of ( O ) 
matters, then it is clear that using an exact diagonaliza- 
tion method for solving a 127 x 127 matrix a number of 
10000 times implies a waste of computer resources. Of all 
the over hundred eigenvalues, only one is needed. Thus, it 
would be more convenient to have a method that targets 
only the desired state with the same accuracy as the exact 
diagonalization method while not having to compute the 
remaining states. In next section we shall compute this 
absortion peak using the density matrix renormalization 
group method |25j, M, pi. 



III. DMRG FOR EXCITONS: RESULTS 




Block A 



FIG. 5: Superblock decomposition of the lattice associ- 
ated to a general tree graph (no cycles allowed). The 
free site denotes the probing site • in the block splitting 
ABC*. 

pof . The advantages of this quantum mechanical DMRG 
(QM-DMRG) over exact diagonalization models in one- 
exciton calculation are the following: 

1. It is possible to reach lattice sizes that are out 
of the reach for exact diagonalization techniques, 
while keeping the same degree of accuracy in the 
determination of energies and wavefunctions. 

2. It is possible to save both CPU time and space for 
the computations in each sampling. This is spe- 
cially interesting when the number of samplings N$ 
is a large number. 



Block A 




We propose to use the DMRG method as explained 
in |2l| for the computation of one-exciton problems with 
off-diagonal disorder. This method has been used to com- 
pute with high accuracy the low lying energy states in 
quantum mechanical problems pflj , not only in one di- 
mensional lattices but also in two and three dimensions 



FIG. 6: Schematic representation of the sweeping process 
in the DMRG process for dendrimers. The free site is at 
a site with coordination number c = 3. 

These properties make the QM-DMRG method spe- 
cially well-suited for excitonic calculations in dimensions 



higher that one when computation means due to disorder 
and connectivity become quite demanding. 

We have tested this proposal in the interisting case of 
extended and compact dendrimers. With this method we 
can increase the number of generations g for each den- 
drimcr family to sizes not reachible with other methods. 

We shall not dwell upon all the details of the QM- 
DMRG formulation that can be found in |2£| ^9|, p0[ , 
but instead we shall stress the main peculiarities of the 
method when applied to extended and compact den- 
drimers. 



A DMRG for Dendrimer Lattices 

For the sake of concreteness, we shall make the dis- 
cussion in terms of dendrimers with a fixed number of 
wedges w = 3 and connectivity c = 3, motivated by the 
Harigaya model, but the method is equally applicable to 
a generic tree graph. Later, we shall present results when 
the number of wedges is different (w 7^ 3). 

In the standard formulation of the DMRG, the whole 
lattice also called universe U or superblock (SB) is split 
into left (Bl) and right (Br) blocks according to the de- 
composition: U — Bl» Br or U = Bl • •Br. The blocks 
describe the degrees of freedom of the system and the 
environment, which can be represented either by Bl or 
Br depending on which stage of the RG-sweeping pro- 
cess we are. In the case of dendrimers, we find more 
appropriate to use the first decomposition based on one 
site, U = Bl* Br. The single lattice site • connects the 
two blocks and serve as a probe to test the reaction of the 
system degrees of freedom to the coupling to the rest of 
the environment. This probing site is a movable site that 
runs all over the lattice during the sweeping process. 

Let us introduce the following notation for the Hamil- 
tonian matrix elements in (p~2h: 



E i=j 
IJij = \j || i -j ||=1 
otherwise; 



(16) 



where i,j are vectors of integer components in a den- 
drimer lattice. We may use a binary notation for labeling 
these sites. 

Although the planar dendrimer lattice extends all over 
the plane, there is a distinctive feature that makes it 
closer to a one single chain. Namely, given any two 
points, there exists a unique path joining them, unlike a 
real two-dimensional lattice like the square lattice. This 
characteristic is reminiscent of the left-right directions 
in one-dimensional chains and it is the basis for trying 
an application of DMRG here |30j . Next we briefly de- 
scribe the several ingredients entering in the application 
of DMRG to dendrimer lattices. 
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FIG. 7: Results for the first absortion peak E a b as in (|l5| ) 
for extended dendrimers versus the generation number g. 
Two numerical methods are compared with the experi- 
mental results of Kopelman et al. DMRG (ours) 
and ED (Harigaya fl22fl ). 

1 Superblock Decomposition of the Dendrimer 

The standard DMRG decomposition of the original lat- 
tice to perform the truncation/renormalization process is 
U = Bf J »»B^ r _ l _ 2 where the subscripts I and N—l — 2 de- 
note the number of sites inside each left and right blocks 
such that their sum - plus two - equals the total num- 
ber of lattice sites N. The index I denotes the iteration 
step of the RG process. This decomposition applies both 
to one-dimensional lattices as well as higher dimensional 
cases p0[ . 

In our case, for a generic dendrimer with a number w of 
wedges and sites with connectivity c, we use a superblock 
decomposition as follows 



U = B 1 {l)B 2 {l)...B c (l)* 



(17) 



where the probing site is joint to the wedge-blocks (see 
Fig||) and it is a movable point that traverses the den- 
drimer lattice as the sweeping process takes places. We 
denote it as a free point in Figj^ and it is a generic point 
of the lattice for a given RG-step I, not always the fo- 
cal point at the core. This decomposition is generic of a 
tree graph with no cycles and the DMRG decomposition 
is applicable with generality to any graph, as shown in 
Fig 
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FIG. 8: DMRG results for the first absortion peak E ah FIG. 9: DMRG results for the absortion peaks E ah (jl|) 
(|l5|) in extended dendrimers versus the generation num- in extended dendrimers (vertical versus the generation 

number g. They corresponds to the ground state and 
two excited states. 

2 Wave-Function Variational Ansatz 



ber g 



Let us assume that the total number of states that we 
want to keep during the truncation/renormalization pro- 
cess is Ne, which includes the GS and Ne — 1 excited 
states. These are the targeted states and we are consid- 
ering the more general case of finding not only the GS 
but also a few of the low lying excited states. Although 
this is not necessary for the one-exciton model that we 
are considering here, however we want to show that the 
method is powerful enough so as to handle more general 
situations. 

The superblock decomposition of the lattice in turn 
induces a decomposition of the wave-function of the tar- 
geted states associated to the blocks and the free site. 
Let us assume extended and compact dendrimers with 
a generic value of w, then the superblock wave-function 
^/(n) at RG-step I and lattice point n is split into the 
following four pieces when the connectivity is c = 3 (17): 



Eafi a a L?(n;l) 
J2 a =i a N E + a L% (n; I) 
J2 a =i a 2N E +aL% (n; I) 

^3N E + 1 



*l(n) 



n 6 Bi(l) 
neB 2 (l) 
n G B 3 (l) 
n = 



(18) 



where denote the free point in Figs.||,|] connecting the 
blocks B c (l),c= 1,2,3; {£°(n; Z)}£f 1( c = 1,2,3 are or- 
thonormal basis of states describing the degrees of free- 
dom in the wedge-blocks, i.e., 



(19) 



and the free unknown coefficients {a a }a=i +1 a * this l- 
stage will be determined later on by means of a diagonal- 
ization/truncation process defining the renormalization. 
They are normalized as 



a r= 



3JV B +1 

E « 



= i 



(20) 



In a generic case of a dendrimer with w wedges, there 
will be a number of blocks c at each stage I of the DMRG, 
and a number of cNe + 1 variational a-parameters. For 
compact dendrimers with w = 3, all ansatzs are like in 
( |l8| ) for every site. However, for extended dendrimers 
with w — 3 we have to distinguish two types of sites, 
those with c = 2 and those with c = 2 (see Fig. ||) In this 
case we need to use another ansatz like in (|18|) for only 



two blocks in (17) when dealing with sites of connectivity 
c = 2. 



3 Superblock Hamiltonians 

The dimensionality of these SB Hamiltonians follows 
from the wavefunction in (|is|), and it is much smaller 
than that of the original Hamiltonian ( |l6| ) . This makes 
their diagonalization something much less demanding 



/V y 

• • i * 

A k 
A. \ 


J-*-* * • * • 

• 


~r\ a 
-A ; \ - 

' A " 


\ A v 

<x / \ 

/ T T 

'A AV 1 



E h (cm ) 



FIG. 10: Pictorical representation of the DMRG wave- 
function for one excited state in an extended dendrimer 
with generation g = 7. The two colors denote positive 
and negative values of the state coefficients and their size 
are proportional to their normalized values. 
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FIG. 11: DMRG results for the first absortion peak 
( |l5| ) in extended dendrimers versus the generation num- 
ber g and with varying number of wedges w. 



4 Truncation of Hilbert Space and Renormalization 



than diagonalizing the whole Hamiltonian. The dimen- 
sion of H$b depends on the number of targeted states. 
In the version of the DMRG we are using here, this di- 
mension is (3Ne + 1) x (3Ne + 1) for both families of 
dendrimers. For example, when the connectivity of a 
given site is c = 3 (where the free site (see Fig.|| ||) lies 
at a given DMRG-step), then the SB Hamiltonian reads 
as follows 
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where H c (l),c — 1,2,3 are the block Hamiltonians for 
each block in the DMRG decomposition of the lattice 
(see Fig.||,||) with dimesions Ne x Ne, namely, 



H 



a,0 



(L°?\H\LP) 



1,2,3; 



(22) 



while H, ltl denotes the value of the original Hamilto- 
nian (|l^) at the free site •/ and the column vectors 
v c (l), c= 1, 2, 3 with dimension Ne x 1 corresponds to the 
interactions between the blocks and the free site. Their 
values depend on the DMRG-step I and constantly up- 
dated through the RG process. The zeroes in © reflect 
the short-range structure of (|l6|). 



The next step is the diagonalization of the superblock 
Hamiltonians in order to obtain the wave functions of 
the Ne targeted states and their energies. The free pa- 
rameters a are then constructed out of the components 
of the targeted wave functions. With this information 
we can perform the projection of the wave functions onto 
the several blocks as well as the renormalization of the 
matrix elements of the superblock Hamiltonian. 

Let us denote the variational parameters by the 
set {aj,a2,a3,«4 ,}i=i> where a*,c = 1,2,3 are Ne- 
dimensional vectors. Assuming that we have three blocks 
(c = 3), then the truncation of the Hilbert space is per- 
formed by the projection of the superblock wave- function 
( p8| ) onto the blocks formed by two out of three possible 
blocks plus the free site. For defineteness, let us assume 
that we are renormalizing blocks B\(l), B 2 (l) at a given 
DMRG-step I, while the remaining block Bs(l) is not di- 
rectly renormalized but only updated through the values 
taken from a previous DMRG-step. Then, the projection 
onto the blocks Bx(l), B 2 (l) is given by 



a = 



{4\ 

4 



(23) 



The projected wave functions in the RHS of ( p3| ) must be 
orthonormalized using a Gram-Schdmit method. They 
will become the new wave functions of the renormalized 
blocks B[,B' 2 . More explicitly, we have 
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where is the orthonormalization matrix 
The new wave functions L'{ (/ + 1) for the renormalizcd 
block B 1 to be constructed in the next RG-step I + 1 are 
given by 

L?(n;/ + 1)= ^E„<aif(n;0 n G B 3 (l) (25) 

After the truncation process we need to renormalize (up- 
date) the different matrix elements of the superblock 
Hamiltonian (|2l| ) . For the blocks B 1 , B 2 we have 

(H X {1) Wl (0\ 
H[(l + 1)=A'U H x (l) v 2 (l)\A' (26) 
\v\(l) vt(l) H. U .J 

/a? ... af-\ 
A' = a'2 1 ... af E (27) 

\< •■■ <v 

These considerations for the renormalization of blocks 
B%(1), -Ba(0 are generic and they equally apply to the 
other pairs of blocks. Moreover, when the free site is 
placed at a lattice site with connectivity c = 2, then the 
renormalization is simplified for then there are only 2 
blocks and the RG process is similar to a one-dimensional 
DMRG step, with a lattice decomposition of the type 
B • B. 



5 Sweeping Dendrimers 

This part of the finite-size DMRG method corresponds 
to moving the free site in Fig. || throughout the lattice 
at each step of the RG process and updating the content 
of the blocks in the DMRG decomposition of the lat- 
tice following the previous renormalization prescriptions. 
The sweeping process is responsible for the convergence 
of the method: the sweeping continues until the energy 
values of the intermediate diagonalizations achieve a de- 
sired prescribed precision. 

During a generic DMRG sweep, the free site starts at 
the focal point and then moves first through the sites in 
wedge w = 1, next the sites of w = 2 are visited and so on 
and so forth with the rest of the wedges until the free site 
returns to the focal point. The sweeping process needs 
a first input called warm-up in order to start up. In the 
case of dendrimers we find useful to do it with a simple 



block renormalization group method (BRG) . This means 
that the dendrimer lattice is decomposed in blocks. The 
blocks in this BRG are formed with 3 sites, one of these 
sites is the father and the other two are children. The first 
step of this BRG starts with a blocking for the outermost 
sites located at the surface of the dendrimer, namely, the 
last and next-to-last generations. During the following 
BRG-steps, the dendrimer lattice is reconstructed from 
the outside to the focal point located at the core. At 
this moment, the warm-up process is completed. In- 
terestingly enough, this form of doing the warm-up in 
the DMRG for dendrimres is similar to one of the ac- 
tual methods of chemical synthesization: the convergent 
method JT|. 

B Numerical Results 

Let us start presenting our numerical results for the 
family of extended dendrimers. Our main purpose is to 
show the capabilities of the DMRG method when deal- 
ing with one-exciton problems in the presence of disorder 
such the case of the Harigaya model that we have taken 
as our model Hamiltonian ( |l2| , [l3"l) . 

We have targeted a number Ne — 3 of states that 
include the ground state plus two excited states. The 
convergence criterion used is to sweep until a precision 
(relative error) of 10~ 10 is achive in all the eigenvalues. 
We find that the degree of convergence with our DMRG 
is very fast. For instance, the number of sweeps needed 
to achieve this accuracy is always lower than 4. The 
number of samples is also Ns — 10000 as in |2^, ||^, 
. The corresponding CPU times are also very small as 
compared with exact diagonalization times. For example, 
for the sampling considered here a typical time for the 
dendrimers with higher generation number is around 900 
s |32j. We start with the number of wedges w = 3. 

In Fig. we present the results for the first absortion 
peak i? a b (^||) in extended dendrimers versus the gen- 
eration number g until a value of 7. The energies are 
measured in cm -1 . Here we are doing two types of com- 
parisons. Firstly, we are comparing our DMRG results 
with the exact diagonalization results (ED) of Harigaya 
up to 5 generations. We clearly see that there is an over- 
lapping of both results and furthermore, we can go easily 
to higher generations. Secondly, we must note in passing 
that there is qualitative agreement between the experi- 
mental results of Kopelman et al. and those of Hari- 
gaya's, although a better quantitative matching would be 
desireable. 

We have been able to reach a generation number of 
g = 13 without using big computer facilities []32| as shown 
in Fig. |[ The vertical lines at each result in this figure 
are error bars due to statistical fluctuations in the sam- 
pling. It is apparent from these results that there is a 
decreasing of the absortion peak as the generation num- 
ber increases, in qualitative agreement with the experi- 



mental results of Kopelman et al. . Although in jlfj 
only extended and compact dendrimers up to g = 5 were 
considered, it has been possible experimentally to reach 
higher generation numbers in a variety of dendrimers of 
different compositions. In our case, we are also inter- 
ested in studying theoretically the behaviour with g of 
physical properties in dendrimers. As it happens, due 
to the exponential growth of the number of sites with g 7 
the dendrimer tends to reach an upper generation limit 
because the surface groups become densely packed. This 
is known as the "de Gennes dense packing". For some 
dendrimers, this limit has been found around g = 10 to 
9 = 12. 

As a result of having this many values we can make a 
polynomial fit to stimate the degree of this slowing down. 
We find that the following scaling law with the generation 
number g 



E c ah = 32000 - C e g a ' cnT 



(28) 

fits well the numerical data with the values C e = 2020 ± 
200cm _1 and the scaling exponent a e = 0.71 ±0.04. The 
continuous line in Fig. ^| gives the fitting in (p8|). 

As we are interested in the performance of the DMRG 
method when dealing with excitons, we also plot two 
byproducts of our numerical calculations. On one hand, 
in Fig^l we present the results for the absortion peaks 
E a b (|l5) corresponding not only to the ground state, but 
also to the lowest two excited states. They are computed 
simultaneously when targeting Ne — 3 states. We see 
that all of these states have decreasing energies as the 
generation number increase. Within the approximations 
of the Harigaya model, this is an indication that these 
lowest lying states are delocalized. On the other hand, 
in Fig. 10 we show a plot of the wavefunction correspond- 
ing to an excited state for a generation number of g = 7. 
This is because the DMRG not only gives energies but 
also computes the corresponding states. In this pictorical 
view of the excited state in the extended dendrimer, each 
circle has a size proportional to the value of the coefficient 
of the normalized wavefunction and the two colors cor- 
responds to positive and negative values. Notice that in 
one of the wedges the state has vanishing contributions. 

So far we have considered the number of wedges fixed 
to w — 3, following the experimental studies in |l6| ] 
and the computations in (22[ M. However, as 
for several types of dendrimers the number of wedges 
achieved ex p or imcnt ally can be different than this num- 
ber |33|, |34|, |35| , we have also applied our DMRG studies 
to deal with extended dendrimers with variable number 
of wedges ranging from io = 2toro = 6as shown in 
Fig. O. In this figure we present the first absortion peaks 
Eah (fl5|) as a function of g and with varying number w. 
It is clearly seen that the effect of increasing w is to lower 
the energy peaks for any generation number (there are 
no level crossings). For a fixed value of w, the behaviour 
of the peaks are all qualitatively similar. 
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FIG. 12: Results for the first absortion peak E ab ( [l5|) us- 
ing DMRG in compact dendrimers versus the generation 
number gpM. 



For the sake of completeness, we have also performed a 
DMRG study in the family of compact dendrimers with a 
fixed number of wedges w = 3 using the Harigaya model 
( |l2| ) , (|ll]) . We plot the first absortion peak E ah ( |l5| ) as a 
function of g up to 8, with the corresponding error bars. 
From this figure it is shown that the energy peaks also de- 
crease with increasing sizes for this family of dendrimers 
and thus the Harigaya model is not able to predict the 
correct experimental data in jl(|. We can make a poly- 
nomial fit to stimate again the degree of this decay. We 
find that the following scaling law 



El b = 32000 - C c g a " cm" 



(29) 

fits very well the numerical data with the values C c — 
3190 ± 120cm _1 and the scaling ex pon ent a c — 0.54 ± 
0.02. The continuous line in Fig. [L2| gives the fitting 
in (|2^). Thus, the only difference that we can capture 
between both families of dendrimers is the fact that the 
scaling exponent a e in the extended dendrimers is bigger 
than a c for the compact ones. 



IV. CONCLUSIONS AND PROSPECTS 

The field of dendrimers has become very active in 
Chemistry during the last decade due to their many pos- 
sible applications explained in the introduction and we 
believe that there are open physical problems and new 
physics so as to make also an entrance in condensed mat- 
ter also. 
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In this work we have been interested in the one-exciton 
processes pertaining to the application of extended den- 
drimers as light-harvesting molecules, i.e., molecular an- n 
tennas. Our main contribution is to introduce the DMRG 
method as a very appropriate computational tool for one- 
exciton approximations with off-diagonal disorder. This 
is because the renormalization/truncation process in a [2 
DMRG calculation is more efficient than the exact diag- 
onalization method if we are only interested in the low- 
lying states of the model Hamiltonian. The gain of the [3 
DMRG method is more pronounced when a large number 
of diagonalizations are needed such the large samplings 
typical of disorder problems with excitons. 

We have shown the potentialities of the DMRG method [4 
in the field of excitons with disorder by computing a vari- 
ety of dendrimer configurations with a variable number of [5 
generations g and wedges w (|l|). In all our computations 
we reproduce the exact results for low generation number jg 
and further, we can increase g without much computa- 
tional effort to numbers that become unreachable for ex- 
act diagonalizations techniques (ED) . Thus we propose to jy 
substitute these ED methods by the DMRG in situations 
like two-dimensional square lattices where the computa- 
tion of one-exciton results with disorder and long-range 
interactions becomes very demanding for ED methods. 
For instance, we have computed lattices with a number 
of 36823 sites in the case of extended dendrimers, just 
to give an example of the sizes that can be reached with 



DM RG. 

One of the open questions m Che ficM of ex- 
tended/compact dendrimers is the explanation of the sce- 
nario in (^) regarding the closing or not of the gap in 
the absortion spectrum using a more elaborate model 
Hamiltonian than the oversimplified Harigaya model 
p2| , p3| |24j ]. It is plausible to think that a truly many- 
body Hamiltonian would be necessary in any attemp to 
address this issue. If we concentrate only in the electronic 
properties of these compounds, one possible candidate to 
describe the physics of multi-excitonic processes is a PPP 
(Pariser-Parr-Pople) Hamiltonian |j, fj], ||. This type of 
Hamiltonian has been used with success in conjugated 
polymers of linear types and we believe that it could also 
be of use for conjugated polymer of dendrimer type. An- 
other simpler candidate for a many-body Hamiltonian 
description of the gap spectrum in extended/compact 
dendrimers is the Frenkel Hamiltonian introduced in (|^) 
that is more manageable for a first try than the PPP 
Hamiltonian. These types of studies are left for future 
work. 



ACKNOWLEDGMENTS 

We would like to thank K. Harigaya for reading a draft 
version of this work. This work has been partially sup- 
ported by the Spanish grant PB98-0685. 



19 

[10 

[11 
[12 
[13 



[14 
[15 
[16[ 

[17[ 
[18 



REFERENCES 

Dendrimer = dendrite + dimcr. From the Greek 
word dendra = tree. Dendrites are the terminations 
of brain cells (neurons), while dimers are bonds (sin- 
glets). 

D. A. Tomalia, P. R. Dvornic, Nature 372, 617 
(1994); R. F. Service, Science267, 458, (1995); J. 
Alper, Chem. Ind. 268, (1991). 

G.R. Newkome, C.N. Moorefield, "Dendrimers" in 
Comprehensive Supramolecular Chemistry; Rein- 
houdt, D. N., Ed.; Elsevier: New York, 1996; Vol. 
10, Chapter 26. 

G.R. Newkome, C.N. Moorefield, C. N.; F. Vogtle, 
Perspectives; VCH: Weinheim, Germany, 1996. 
O.A. Matthews, A.N. Shipway, J.F. Stoddart, Prog. 
Polym. Sci. 1-56, (1998). 

D. Baeriswyl, D.K. Campbell and S. Mazumdar, 
in Conjugated Conducting Polymers, edited by H. 
Kiess (Springer- Verlag, Heidelberg, 1992), pp 7-133. 
S. Pleutin, E. Jeckelmann, MA. Martin-Delgado, 
G. Sierra, "The Recurrent Variational Approach 
applied to the Electronic Structure of Conjugated 
Polymer s", Progress in The oretical Chemistry and 
Physics; |cond-mat/9908062| . 

MA. Martin-Delgado, G. Sierra, S. Pleutin, E. Jeck- 
elmann, " Matrix Product Approach to Conj ugated 



Polymers", Phys. Rev. B 61, 1841 (2000); cond- 



mat/9908066 



C. J. Suckling, J. Chem. Soc, Chem. Commun., 661, 
(1982). 

J. A. Hyatt, J. Org. Chem. 43, 1808-1811, (1978); 
F. Vogtle, E. Weber, Angew. Chem. Int. Ed. Engl. 
13, 814-815, (1974). 

Y. Murakami, A. Nakano, K. Akiyoshi, K. Fukuya, 
J. Chem. Soc, Perkin Trans 1, 2800-2808 (1981). 
E. Buhleier, W. Wehner, F. Vogtle, Synthesis, 155- 
158, (1978). 

D. A. Tomalia et al., Polym. J. 17, 117-132, (1985); 
D. A. Tomalia et al., Macromolecules 19, 2466-2468, 
(1986). For an extensive review, see: D. A. Tomalia, 
A. M. Naylor, W. A. Goddard III, Angew. Chem. 
Int. Ed. Engl. 29, 138-175, (1990). 

R. G. Denkewalter, J. Kolc, W. J. Lukasavage, US 
Pat. 4, 289, (1981). 

C.J. Hawker, J.M. Frechet, J. Am. Chem. Soc. 112, 
7638-7647, (1990). 

R. Kopelman, M. Shortreed, Z-Y. Shi, Z. Xu, J.S. 
Moore, A. Bar-Haim, J. Klafter, "Spectroscopy Evi- 
dence for Excitonic Localization in Fractal Antenna 
Supermolecules" , Phys. Rev. Lett. 78 1239, (1997). 
J.F.G.A. Jansen, E.M.M. de Brabander-van den 
Berg, EW. Meijer, Science 266, 1226-1229 (1994). 
C. Devadoss, P. Bharathi, and J. S. Moore, J. Am. 
Chem. Sec. 118, 9635 (1996); M. Shortreed, Z-Y. 
Shi, and R. Kopelman, Mole. Cryst. Liq. Cryst. 28, 
95 (1996); M. R. Shortreed, S. F. Swallen, Z-Y. Shi, 



12 



W. Tan, Z. Xec, C. Devadoss, J. S. Moore, and R. 
Kopelman, J. Phys. Chem. B. 101, 6318 (1997). 
[19] A. Bar-Haim, J. Klafter, and R. Kopelman, J. Am. 
Chem. Soc. 119, 6197 (1997); A. Bar-Haim and J. 
Klafter, J. Phys. Chem. B 102, 1662 (1998); A. Bar- 
Haim and J. Klafter, J. of Luminescence 76-77, 197 

(1998) ; A. Bar-Haim and J. Klafter, J. Chem. Phys. 
109, 5187 (1998). 

[20] S. Raychaudhuri, Y. Shapir, V. Chernyak, S. 
Mukamel, "Excitonic Funneling in Extended Dcn- 
drimers with Non-Linear and Random Potentials" , 
Phys. Rev. Lett. 85, 282 (2000). 

[21] The excitonic operators 6j are in fact hard-core 
bosons. The local states |vac)j, |exc)i are in one-to- 
one correspondence with spin states |— )j, and 
likewise, the operators bi, b\ and b\bi can be rep- 
resented by Pauli matrices cr~, of and (of + l)/2, 
respectively. Thus, the Hamiltonian (||) becomes the 
XY model with transverse couplings given by Jij 
and external magnetic field given by Ei. 

[22] K. Harigaya, "Coupled exciton model with off- 
diagonal disorder for optical excitations in extended 
dendrimers" Phys. Chem. Chem. Phys. 1, 1687 

(1999) . 

[23] K. Harigaya, "Optical excitations in diphenylacety- 
lene based dendrimers studied by a coupled exciton 
model with off-diagonal disorder " Int. J. Mod. Phys. 
B 13 2531, (1999). 

[24] K. Harigaya, "Optical Excitations in Fractal An- 
tenna Supermolecules: Extended Dendrimers" in 



4th Int. Symposium on Advanced Physical Fields: 

Quantum Phenomena in Advanced Materials at 

High Magnetic Fields. 
[25] S.R. White, Phys. Rev. Lett. 69, 2863 (1992). 
[26] S.R. White, Phys. Rev. B 48, 10345 (1993). 
[27] S.R. White and R. Noack, Phys. Rev. Lett. 68, 3487, 

(1992). 

[28] M.A. Martin-Delgado, G. Sierra and R.M. Noack, J. 

of Phys. A: Math, and Gen. 32, 6079 (1999). 
[29] M.A. Martin-Delgado and G. Sierra, Phys. Rev. 

Lett. 83, 1514 (1999). 
[30] M.A. Martin-Delgado, J. Rodriguez-Laguna, G.Sierra; 

Single-Block Renormalization Group: Quantum 



Mechanical Problems, |cond-mat/0009474| (2000). 

[31] The error bars for g = 12, 13 are bigger in these 
generations just because the sampling number is re- 
duced, namely, N$ = 100. 

[32] We have used a personal computer with a Pentium 
III at 450 MHZ. 

[33] J. M. J. Predict, J. Am. Chem. Soc. 112, 7638-7647 
(1990). 

[34] R. Roy, D. Zanini, S. J. Meunier, A. Romanowska, 
J. Chem. Soc, Chem. Commun. 1869-1872 (1993). 

[35] An example of a 4-lobule dendrimer is the Porphyrin- 
cored third generation dendrimer synthesized in P. 
J. Dandliker et al. Angew. Chem. Int. Ed. Engl. 33, 
1739-1742 (1994); H.-F. Chow et al., Chem. Eur. J. 
2, 1085-1091 (1996). 



